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Abstract 

When Brans-Dicke Theory is formulated in terms of the Jordan 
scalar field (p, dark energy is related to the mass of this field. We in- 
vestigate the solution which is relevant for the late universe. We show 
that if <p is taken to be a complex scalar field then an exact solution 
of the vacuum equations shows that Friedmann equation possesses a 
term, proportional to the inverse sixth power of the scale factor, as 
well as a constant term. Possible interpretations and phenomenologi- 
cal implications of this result are discussed. 
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1 Introduction 



Jordan in 1959 proposed a scalar field 99 and replaced the gravitational con- 
stant in the Einstein- Hilbert action by a term proportional to ip^ [1, 2]. Brans 
and Dicke redefined the scalar field so that □(/> = as a vacuum solution [3] . 
Thus, in Jordan-Brans-Dicke theories (JBD), the counterpart of the gravi- 
tational coupling term, 1/1QtiGm, is replaced by or Lp'^/Su, which may be 
a function of space and/or time. For an isotropic homogeneous cosmology, 
which evolves in time, scalar field is a function of only time. JBD gravity 
has been used extensively to develop dark energy models [4-9] which usually 
involve a scalar potential which is adjusted to fit observed cosmology. In con- 
trast; we will use a simple model where the JBD field is complex and besides 
the kinetic term, it contains only a standard mass term for the Jordan field 
cp. This is in contrast to most cosmological models based on Jordan-Brans- 
Dicke theory where an arbitrary potential V{(fi) is assumed and then the 
potential is adjusted to give the desired solution. One possible motivation 
for considering a phase for the JBD field is that it could affect the strong 
QCD parity violating phase x [10, 11] or the CP violating phase of the quark 
mass matrix [12, 13] or neutrino mixing [14]. These may be related to dark 
matter. While in Section II, an exact cosmological solution is displayed and 
its stability analysis is made for the vacuum case, in section III, it is inter- 
preted as the effect of the proposed complex scalar field to the evolution of 
the late time universe. We also give a brief discussion of possible relationship 
of the phase of the complex scalar field. We will show that the dark energy 
component, in addition to the constant term makes a contribution to the 
Friedmann equation, which is proportional to the inverse sixth power of the 
scale factor. 

We use a metric signature (-1 ) and the Jordan formalism. The 

lagrangian densities of the Jordan and Brans-Dicke language are related by, 

= -^R + \g^''d^ipd,ip = £jBD. (1) 
We add a mass term and a phase x, 

</? = </7i i</?2 = Iv'l e'^ (2) 
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so that the lagrangian density with a complex massive scalar field becomes, 

(3) 



which can also be expressed as 
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The action is defined by 



-I 



ct^x\/—g£ + 5" 



(5) 



When this action is varied with respect to the metric and the complex scalar 
field, the equations of motion, for a Friedmann- Robertson- Walker metric and 
perfect fluid energy-momentum tensor — diag (p, —p, —p, —p) , reduce to 
the following: 
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(f = Q. 
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Equation 8, being complex, is equivalent to two real equations. H and F\-\-iF2 
are respectively defined as the fractional rate of changes of the scale size and 
the JBD scalar field </?. 



a iD 
H=- , Fi + iF2 = - 
a if 



(9) 
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where Fi and F2 are real, and defined as, 



^2 = — = X- 11 

For spatially fiat (k=0) universe, we will obtain solutions which give H, Fi, 
F2 as a function of the scale size a. After rewriting equations 6, 7, and 8 in 
terms of H, Fi, F2, and their derivatives with respect to the scale size of the 
universe, a, we get the following equations: 

3//2 - 2u{Fl + Fi) + mF^ - 2um^ = (12) 



Si/^ + (2a; + 4)Fi^ + 2a;F2 + AHFi + 2aii'(F; + H') - 2um^ = ^^^^ (13) 



-6H^ + 2uF^ - 2uFi + 6uHFi + 2aojHF[ - 3aHH' + 2um^ = (14) 

(4a;Fi + 6a;//) F2 + 2ujaHF^ = (15) 
where prime denotes derivative with respect to a. 



2 The exact cosmological solutions for dust 
dominated and the vacuum case 

Exact cosmological solutions of JBD theory can be obtained by analyzing 
symmetries of the field equations [15]. Perturbative solutions were found in 
[16, 17] . We propose an ansatz for finding the exact solutions for the late time 
universe. We will show that an exact solution of the cosmological equations 
can be obtained by an ansatz for the matter dominated case and vacuum 
case. 
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By integrating equation 17 with the ansatz yields 



a 
ao 



(18) 



Since cu > 10^ [18-21] \ip\ varies very slowly as the universe evolves so ip is 
approximately constant during matter dominated era. Variation of gravita- 
tional constant which is proportional to Fi in equation 10 is also small due 
to its proportionality to K 

By using the ansatz, equation 16, equation 15 gives the solution 
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since cu > 10^ for all practical purposes a — —3. For dust solution of complex 
JBD model, pm = Po j and pu = 0, when equation 16 and 19 are placed 
into equation 12, equations 13, 14 are satisfied. Equation 21 is derived as, 

2a' 



(3c^ + 4)(26^ + 3) 
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2p 



•0 



+ n 



(21) 



Several interesting features emerge for the vacuum case, pM = Pm = 0. 
Then contains only a constant term and a 1/a® term, respectively. Since 
F2 is directly related to the phase of the scalar field in equation 11, the 
term is identified with the effect of the phase to the expansion of the universe, 

H'{a) = , M^ + ^f .2 p2^^^y (22) 
^ ' (3w + 4)(2a; + 3)^ " ^ ' 

When we restrict the model to the F2 = case, F = Fi, complex scalar 
field turns into a real scalar field and previously studied JBD equations are 
obtained [17]. To test the stabihty of our vacuum solution; we set 

a = a(l + erj) (23) 
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</? = (p{l + eijj). 



(24) 



Inserting these variables into Eq. 12-15, and after neglecting the higher order 
terms in e, we obtain four homogeneous differential equations for the three 
functions r],ipR, and ipj. We set 

V = Voe^' (25) 
= ^me^' (26) 

ipi = ipioe^' (27) 

and obtain four linear homogeneous equations for the three unknowns rjo, ipRo, 
and ■0/0- This homogeneous system has a 4x3 matrix of coefficients. The 
condition that a nontrivial solution exists is that the rank of the matrix of 
coefficients is at most two. All 3x3 subdeterminants must be zero to obtain 
a nontrivial solution for /3. 

We thus obtain four equations for the four 3x3 subdeterminants, and 
arrange them in powers of (3. One equation is cubic in (3 whereas the other 
three are quadratic. 

AnP^ + ^12^' + Ais/3 + >li4 = (28) 
+ + ^4 = (29) 
^32^' + AssP + A34 = (30) 

A42/3' + A43/3 + A44 = (31) 

The determinant of the matrix of coefficients (dct A) is nonzero, so we 
can conclude that there is no solution for (3 and it means that the solution 
is stable. 

With the addition of the cosmological constant to cold dark matter(CDM) 
model, the resulting ACDM model gives a better fit [22,23]. Five- Year 
Wilkinson Microwave Anisotropy Probe (WMAP) temperature and polariza- 
tion observations [24] which include data from Baryon Acoustic Oscillations 
in the galaxy and Type la supernova luminosity/time dilation measurements 
[26] are used in the fitting process. 

The constant term in equation 21 plays the role of the cosmological con- 
stant and we will investigate the extra term from the phenomenological point 
of view. 
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3 Phenomenology 



For standard cosmology in "matter dominated" era, 

^ = Qa + Qm(-)' (32) 
rzg a 

where Q/^ is the fraction of vacuum energy and the matter fraction, Qm 
is interpreted as Qm = ^vm + ^dm where Qvm is the fraction of visible 
matter and Qdm is the fraction of dark matter. For the complex JBD model, 
Friedmann equation becomes, 

^)^a..a.(^)%a.(^)'. (33) 

We make the variable transformation and define 

u = VirJ^y + ^ (34) 
\aoJ 



c=,/|Si.-||l (35) 
^ = sgn(n^-^J (36) 
SO that Eq.33 is put in differential form 

3^Hodt = f" (37) 

To test the viability of equation 33 predicted by complex JBD model, we 

have to compare the standard fit to union data [26] with a fit using equation 
33 and Hq=71 km/sec/Mpc. With the constraint Q/^ + Qm = 1, the standard 
model fit has one free parameter which can be chosen as Q\. 

The latest Type la supernova (SNIa) data sets in Table C2 are taken by 
the Supernova Cosmology Project Group [26]. 414 SNIa supernova magnitude- 
redshift observations are compiled with "Union" and after selection cuts, it 
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reduces to 307 Sne. Luminosity distance, in Priedmann-Robertson- Walker 
Cosmology, is defined in [27] 

dl = J dz'[j2m + zT'^""'^ - + z'f]"^ 

where Kg = ~ ^■ 

Distance modulus can be written in terms of luminosity distance 

;, = m-M = 51og(-^) + 25, (38) 
Mpc 

where m is the apparent magnitude and M is the absolute magnitude. Three 
different cases for x are analyzed. Three fits we would like to present are; 

1°) X = —1, Qa=0 standard cosmology fit with Qa + Qm = 1 gives, 

(1-)^ 0.745 + 0.255 (^)' (39) 

with -y^ jd.o.j . = 1.45. The scale factor and lifetime of universe are related 

by, ^ 

^ [cosh{3HoV^t) - 1) (40) 



2°) X = +1 fit to complex JBD model with Qm — and Qa + Qa — 1 
gives, 

(|-)^ 0.938 + 0.062 (:^)' (41) 

with /d.o.f. = 1.43. The scale factor and lifetime of universe are related 

by, 

-) = \l^smh{3HoV^t). (42) 

Oo/ V "A 



3°) X = , fit to complex JBD model with matter, + CIm + ^^A = 1 



q2 

and flA = jn^ so that equation 33 becomes 



/OA + vr^A(-)' (43) 

tin a 
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Figure 1: Magnitude vs. Redshift Graph 

gives, 

(^)^ 0.7.0 + 0.200 (^)%0.010(^)' (44, 

with x^/(i.o./. = 1.43. This condition makes the r.h.s. of Friedmann equation 
perfect square and gives the best fit to supernovae union data. The scale 
factor and hfetime of universe are related by, 

(^)' = if^'"p^^^°^'^"'^ ^^^^ 

The fit for equation 41 is interesting but unacceptable since the model 
estimated lifetime is approximately 10 Gyrs, small compared with the ob- 
servations [25]. In fig.l, standard cosmology fit using Eq.(39) and complex 
JBD fit Eq.(44) to Union data sets [26] arc shown. Wc conclude that the 
extra term, from the change of the JBD field's phasc,which must be small, 
indicates a shghtly better fit. The difference between two models will be 
seen for high redshift observations. Similar to the addition of cosmological 
constant in ACDM, cosmological reason of obtaining the phase term in JBD 
theory can be investigated. 
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4 Conclusion and Acknowledgements 



Complex JBD model which in addition to the constant term makes a con- 
tribution of ^ term to the Friedmann Equation fits the Supernovae data 
accurately. It is clear from 22, which is valid for pM = Pm = 0, that this 
term is a natural component of dark energy. Actually we conclude that the 
complex JBD model explains the evolution of the universe with a slightly 
better fit. However; the lifetime of the universe was found approximately 
12 Gyrs, smaller than the experimental age of the universe, 13,8 Gyrs [25]. 
Moreover; the complex scalar field is responsible in large scale for the ex- 
pansion behavior of universe and its phase behaves as a phenomenological 
extreme density term. The more the density due to complex phase becomes, 
a smaller model-based age is determined. In standard cosmology, the ^ 
term would be obtained if the kinetic term of a scalar field dominates the 
energy-momentum tensor. In our model, this is also how it mathematically 
arises. However, the physical interpretation is that it IS; dbS db typical feature 
of CJBD theory, a component of dark energy and its presence may be de- 
termined by more accurate measurements of the Supernova data near z ~ 3. 
This is far from the radiation dominated age, z ^ 1100. This research is in 
part supported by The Turkish Academy of Sciences, TUBA. 
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